Fulde-Ferrel-Larkin-Ovchinnikov Inhomogeneous Superconducting State and Phase 
Transitions Induced by Spin Accumulation in a Ferromagnet-d^ ,^-Wave 
Superconductor-Ferromagnet Tunnel Junction 
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Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) inhomogeneous superconducting (SC) state, first- and 
second-order phase transitions, and quantum criticality induced by spin accumulation in a 
ferromagnet-d :E 2_j / 2-wave superconductor-ferromagnet tunnel junction are theoretically predicted. 
A complex phase diagram in the temperature-bias voltage plane is determined. It is found that the 
phase transitions from the homogeneous BCS state to the inhomogeneous FFLO state, and from 
the FFLO state with the momentum q's azimuthal angle # q = to that with # q = 7r/4, are of the 
first-order; while the transitions from all SC states to the normal state at critical voltages are of the 
second-order. A Lifshitz point, a bicritical point and a quantum critical point are identified. 

PACS numbers: 74.50.+r, 75.47.De, 73.40.Rw 



Introduction. — Fortyyears ago Fulde and FerrelQ], 
Larkin and Ovchinnikov 2J (FFLO) pointed out indepen- 
dently that owing to the presence of the spin-exchange 
field, the superconducting (SC) order parameter in a fer- 
romagnetic superconductor can be spatially modulated 
in real space, leading to that the SC state inside these 
materials is inhomogeneous. Such a spatially modulated 
SC state has long been detected in a number of systems 
under various circumstances (e.g. 0,0). Nonetheless, 
a clear and direct confirmation of the existence of the 
FFLO state has not yet been reported sofar, because a 
successful observation of the FFLO state requires very 
clean type II superconductor with the characteristic pa- 
rameter k y§> 1 , the condition is not easily satisfied in 
conventional superconductors. Recently, the possibility 
of the FFLO states in ultracold atomic Fermi gases has 
been actively discussed (see e.g. Refs. and a new 

phase observed in CeCoIn 5 is also found to be the FFLO 
state with a spatially modulated superconducting order 
parameterQ. On the other hand, for a magnetic sand- 
wiched heterostructure with a normal metal or super- 
conductor as a spacer, when the spin-polarized electrons 
driven by an external bias voltage from a ferromagnetic 
film enter into the spacer faster than the spin-polarization 
can diffuse away from the interface, there must be some 
nonequilibrium spin densities built up near the interfaces 
of the heterostructures that depend on the relative orien- 
tations of magnetic moments, leading to the occurrence 
of spin accumulation. The spin accumulation is an im- 
portant effect in the emer ging fie ld of spintronics (see, e.g. 
Refs. H H [H EH E El 111 for review) , which plays 
a key role in the spin Hall effect ^3 as well as in various 
spintronic devices such as spin transistors, spin valves, 
spin diodes, spin field-effect transistor, and so on. Quite 
recently, the spin accumulation has been experimen tally 
observed in a n-type GaAs by means of Kerr rotation 
Another interesting observation is that the spin accumu- 
lation could result in the suppression of superconductiv- 



ity, that happens in ferromagnet /superconductor (F/S) 
heterostructures (e.g. Refs. [H GJ, Hi |H HJ IH ) . 
The cause is that the injected spin-polarized electrons 
into the central superconductor give rise to a spin im- 
balance, leading to a spin density accumulated near the 
interfaces. Such a spin density is equivalent to a small 
magnetic field that plays a role as a pair-breaking field, 
thereby enabling the superconductivity to be suppressed. 
As the spin accumulation in the central superconductor 
could generate an equivalent magnetic field, one has very 
reasons to expect that the FFLO spatially modulated SC 
state should appear in such a system. 

There has been a recent study|2J] showing that in a 
certain range of bias voltage the FFLO inhomogeneous 
SC state can indeed occur in the central superconductor 
of a ferromagnet/s-wave superconductor double tunnel 
junction in case of antiparallel alignments of the magne- 
tizations of both ferromagnets. It has been found that 
the phase transition from the homogeneous s-wave BCS 
state to the inhomogeneous FFLO state is of first-order, 
that gives rise to the oscillating behaviors of the tunnel 
magnetoresistance (TMR) and conductance. These os- 
cillating features can be used as an alternative way to 
examine the existence of the FFLO state. In this paper, 
we shall investigate the FFLO state induced by the spin 
accumulation in a central d^^^-wave superconductor 
of the ferromagnet/d-wave superconductor double tun- 
nel junctions. As the SC gap in a d-wave superconductor 
is anisotropic, it would induce more exotic behaviors than 
in a s-wave superconductor. It is shown that the phase 
transitions from the homogeneous d x 2_ y 2-w&ve BCS state 
to the inhomogeneous FFLO state, and from the FFLO 
state with the azimuthal angle of the momentum # q = 
to that with # q = ir/4, are of the first-order; while the 
transitions from all SC states to the normal state at crit- 
ical voltages are of the second-order. A complex phase 
diagram in the temperature-bias voltage plane is deter- 
mined, and a Lifshitz point, a bicritical point and a quan- 
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FIG. 1: (Color online) The schematic illustration of the F/d- 
SC/F double tunnel junction. The lower depiction indicates 
the d x 2_ y 2 -wave gap symmetry in the central superconductor, 
with 8q the azimuthal angle of the momentum q. 



turn critical point are identified. 

Model and formalism. — We start with a symmetric 
ferromagnet jd x i_ y i -wave superconductor/ferromagnet 
(F/d-SC/F) double tunnel junction with the left and 
right ferromagnetic (FM) electrodes applied by bias 
voltages V/2 and —V/2, respectively. The two identical 
FM electrodes with antiparallel alignment of both 
magnetizations are separated from the central d x 2_ y 2- 
wave superconductor by two insulating thin films. The 
d x 2_ y 2-w&ve superconductor is presumably described by 
the BCS framework. Suppose that the energy relaxation 
time of quasiparticles in the superconductor is shorter 
than the time between two successive tunneling events, 
while the latter is shorter than the spin relaxation 
time, which ensures that the electrons tunneling into 
the superconducting spacer could comply the Fermi 
distribution, and meanwhile keep their spin directions. 
For simplicity, the Andreev reflection effect will be 
reasonably ignored, as the resistance of this tunnel 
junction with insulating thin films is presumed to be 
greater than that of a conventional metallic contact. 
The schematic layout of the system is shown in Fig. 1. 

In terms of the standard tunneling Hamiltonian and 
within the linear response theory, the tunne ling current 
through the zth junction can be obtained by[ 
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In these above equations, fo(z) is the Fermi distribu- 
tion function of thermal equilibrium, and f^ a {a =],[) 
is the noncquilibrium distribution function of quasiparti- 
cles. The excitation energy is given by E^ = v^k + ^k 
for the central d x 2_ y 2-w&ve superconductor. The d x 2_ y 2- 
wave SC gap parameter is defined as usual by Ak = 
A cos(26>k), where 6*k = tan -1 (k y /k x ) is the azimuthal 
angle of the momentum k, and £k is the single-particle 
energy measured from the Fermi level. The coher- 
ence factors Uk and «k are determined by w k = i(l + 



At). 5* and 



A*), and vl = Kl-&/v^k-r-k 
Q denote the spin density and the quasiparticle charge 
density, respectively. 

FFLO state in a d x 2_ y 2-wave superconductor. 
As mentioned above, in the case of antiparallel align- 
ment of magnetizations, the nonequilibrium spin accu- 
mulation near the interfaces of the tunnel junction can 
appear [ItI I2I ]. thereby giving rise to a chemical po- 
tential shift that plays essentially the same role as the 
spin-exchange field explored by FFLO in their seminal 
articles^ U Since the tunnel junction under interest 
is symmetrical, we may presume the chemical poten- 
tial shift (5rt| = -<5aH = H~|- Generally, the 
FFLO state involves the finite momentum q pairing of 
electrons [2^. In the following, we will consider the FFLO 
state just in the sense of the Cooper pairing with nonzero 
momentum q. For a d x 2_ y 2-w&ve gap symmetry with fi- 
nite momentum q pairing, the energy of quasiparticles 
in the presence of the chemical potential shift induced 
by the nonequilibrium spin accumulation has the form of 

E k a = y/& + A k + CO s(0k- q )- Sfi], where v F 

is the Fermi velocity. The amplitude A is determined by 
the gap equation 
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where T is the tunneling matrix element, i = 1, 2, Di„ is 
the subband density of states (DOS) in the iih FM elec- 
trode, a = ±1 for spin up and down, respectively, rji = 1, 
T]2 = — 1, S, Q, N and N are given by 



where Vbcs is the BCS-type pair interaction, N(0) is 
the DOS in the normal state at the Fermi level, and w c 
is the cutoff energy. In the absence of spin-flip scatter- 
ing, the spin up and down tunneling currents should be 
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conserved, i.e. I\ a — 12a, leading to 

S = PN, (7) 
0=0, (8) 

where P = \D^ ~ Dii\/\Dii +Dii \ is the spin polarization 
of the FM electrodes. Eq. JSJ tells that there is no charge 
accumulation in the system. Following the treatments in 

Refs.mHil, 

we consider the solution of form 

fi ar = ME k -*8n). (9) 

It is applicable when the thickness of the central super- 
conductor is much smaller than the spin diffusion length, 
and the spin relaxation time is sufficiently long|2£j. 
These above equations should be solved self-consistently. 

As discussed by Abrikosovj^, for an inhomogeneous 
superconducting state, the nonzero solution for A implies 
only the local minimum of the free energy. For the multi- 
ple solutions for different q, only the value of q (thus # q ) 
that gives the lowest free energy of the system is kept. 
The free energy of the present system bears the following 
formic 

(10) 

where F$ and F N represent the free energy of the SC 
state and the normal state, respectively. By comparing 
the free energies of different phases, the phase diagram 
of the central superconductor can be obtained. 

Phase diagram. — Figure 2 presents the phase diagram 
of the central d x 2_ y 2-wa,ve superconductor in the tem- 
perature (T/T c )-bias voltage (eV/2A ) plane for P = 
0.4, where T c is the SC transition temperature, and 
Ao = A(T = 0). Along V c line, the superconductiv- 
ity is completely suppressed. It turns out that the V c 
line is the phase boundary that separates the SC states 
from the normal state, and consequently, the phase tran- 
sition along the V c line is of second-order. Below the 
V c line, there are three phases in four regions: (i) the 
d-wave BCS phase, which is in a homogeneous d x i_ y i- 
wave BCS superconducting state with zero momentum 
pairing; (ii) the region encompassed by CiB 1 B 2 Ci, which 
is in an inhomogeneous FFLO state with finite momen- 
tum pairing at # q = 0; (iii) the region encompassed by 
C1B2B3C2LPC1, which is in the inhomogeneous FFLO 
state with finite momentum pairing at # q = 7r/4; (iv) 
the small region encompassed by C2B3VQCPC2, which 
is in the FFLO state at # q = 0, the same phase as in 
(ii). The phase boundaries, LPC1B1, CiB 2 , and C2B3, 
separate these four regions, along which the phase tran- 
sitions between the homogeneous d-wave BCS state and 
the inhomogeneous FFLO states with 6> q = and 7r/4 are 
of the first-order. This statement is clearly confirmed in 
Fig. 3, where the chemical potential shift induced by the 




FIG. 2: (Color online) The phase diagram of the central 
d x 2_ y 2-wave superconductor in the temperature-bias volt- 
age plane for the antiparallel alignment of magnetizations at 
P = 0.4. Four phases are observed: the homogeneous d x 2_ y 2- 
wave BCS state with zero momentum pairing, the FFLO 
states with finite momentum pairing q for 6* q = and 7r/4, 
and the normal state. A Lifshitz point and a bicritical point 
are identified, and a quantum phase transition is sepcified at 
Vqcp- 



spin accumulation versus the bias voltage shows two dis- 
.cojatinvjities for 6* q = and 7r/4 at specific voltages. The 
inset of Fig. 3 presents the gap amplitude A/Ao ver- 
sus the bias voltage eV/2Ao, which also shows that the 
transitions from the d-wave BCS state to the FFLO state 
with 6* q = 0, and then to the FFLO state with q = 7r/4 
are discontinuous. These discontinuities characterize the 
first-order phase transitions between those different SC 
states. With increasing the bias voltage, the chemical 
potential shift is increasing in trend and coincides with 
that in normal state, while the gap amplitude is decreas- 
ing and vanishes at V c , as shown in Fig. 3. Our calcula- 
tions demonstrate that no other stable FFLO states with 
8q other than and 7r/4 can be found. In another word, 
the FFLO phases with # q = and 7r/4 have the lowest 
free energies than those with any other values of # q . 

In the phase diagram, there are three characteristic 
points along V c line: LP, C2 and Vqcp- At the point 
LP, a disordered phase (normal state) , an ordered phase 
(d-wave BCS state) and a spatially modulated phase 
(FFLO state with 6* q = 7r/4) meet, showing that the 
point LP is a Lifshitz point 29]. Although the Lifshitz 
point was originally proposed for a metamagnet, we now 
have another example in a different system. At the 
point C2, two second-order phase transition lines meet 
with a first-order phase transition line, suggesting that 
C2 is a bicritical point. At the point VqcPj satisfying 
PsVqcp/2Ao — 1.03, a quantum phase transition (QPT) 
from the ordered FFLO state with # q = to the disor- 
dered normal state happens at T = 0, indicating that 
Vqcp is a quantum critical point. Such a QPT is of 
second-order. 

At the point Ci, three first-order phase transition lines 
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FIG. 3: (Color online) The bias dependence of the chemical 
potential shift induced by the spin accumulation at P — 0.4, 
T/T c — 0.2, where two discontinuities are observed. Inset: 
the bias voltage dependence of the gap amplitude (A/Ao). 



FIG. 4: (Color online) The bias voltage dependence of the 
tunnel conductance at P — 0.4, T/T c — 0.2, where N denotes 
the normal state. 



meet, implying that the three SC phases can coexist. 
Thus, Ci is a triple transition point. At the points Bi, 
B2, B3, the first-order phase transitions may appear at 
T = 0, namely, at Bi the homogeneous d-wave BCS 
state coexists with the inhomogeneous FFLO state with 
9q = 0; and at B 2 , B 3 , the FFLO states with # q = and 
6q = tt/4 can coexist in the ground state. These first- 
order phase transitions at zero temperature are driven by 
external voltages, not by thermal fluctuations. 

Tunnel conductance and magnetoresistance. — The 
tunnel currents are given by I A = Io(l — P 2 )N for 
the antiparallel (A) alignment of magnetizations, and 
I p = IqN for the parallel (P) alignment of magneti- 



zations, where Iq = 2ire 



T 



(At + Ai). 



The tun- 

dI A > p /dV, 



nel conductance can be obtained by G A ' P 
and the tunnel magnetoresistance (TMR) is defined by 
TMR = — 1. Fig. 4 gives the tunnel conduc- 
tance G A ' P /Gjf versus the bias voltage (eV/2A ) at 
T/T c = 0.2 for P = 0.4, where G^ is the tunnel conduc- 
tance in the normal state, and the curves of the conduc- 
tance corresponding to different SC phases are indicated. 
It can be seen that the conductance G A in the A config- 
uration shows exotic behaviors with jumps and kinks, in 
comparison to the conductance G P in the P configura- 
tion, which reflects the features of the superconducting 
density of states. This behavior is also different from the 
s- wave superconductor explored in Ref. ji^l . Fig. 5 shows 
the bias voltage dependence of the tunnel TMR. One may 
observe that the TMR displays complex behaviors with 
cusps and kinks, being resulted from the first-order phase 
transitions between different SC phases. The behavior of 
TMR is also obviously different from the s-wave super- 
conductor discussed in Ref. ■ The results of the tunnel 
conductance and TMR show that the anisotropy of the 
SC gap parameter indeed has an essential effect on the 
spin-dependent transport of electrons in F/SC/F double 
tunnel junctions. 

Summary. — In summary, we have theoretically in- 
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FIG. 5: The bias voltage dependence of the tunnel magne- 
toresistance at P = 0.4, T/T c = 0.2, where N denotes the 
normal state. 



vestigated the effect of spin accumulation on the spin- 
dependent transport in F/d-SC/F double tunnel junc- 
tions. The FFLO inhomogeneous SC state, first- and 
second-order phase transitions, and quantum criticality 
induced by the spin accumulation in the central d x i_ y i- 
wave superconductor are predicted. A complex phase 
diagram in the temperature-bias voltage plane is deter- 
mined. It is shown that the phase transitions from the 
homogeneous BCS state with zero momentum pairing to 
the inhomogeneous FFLO state with finite momentum 
pairing, and from the FFLO state with 6> q = to that 
with # q = 7r/4, are of the first-order; while the transitions 
from all SC states to the normal state at critical voltages 
V c are of the second-order. A Lifshitz point, a bicrit- 
ical point and a quantum critical point are identified. 
The tunnel conductance and TMR are also obtained. 
It is found that the effect of the spin accumulation on 
the spin-dependent transport in F/d-SC/F double tun- 
nel junctions is quite different from that in F/s-SC/F 
tunnel junctions. 

This work is supported in part by the NSFC (Grant 
Nos. 90403036, 20490210, 10247002). 
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